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Inspired by the multicanonical approach to simulations of first-order phase transitions we propose for g-state 
Potts models a combination of cluster updates with reweighting of the bond configurations in the Fortuin- 
Kastelein-Swendsen-Wang representation of this model. Numerical tests for the two-dimensional models with 
q = 7, 10 and 20 show that the autocorrelation times of this algorithm grow with the system size V as t <x V a , 
where the exponent takes the optimal random walk value of a « 1. 
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1. INTRODUCTION 

Monte Carlo simulations of first-order phase 
transitions [1] in the canonical ensemble are 
severely hampered by extremely large autocorre- 
lation times r oc exp(2crL D ~ 1 ) where a is the 
(reduced) interface tension between the coexist- 
ing phases and L D ~ X is the cross-section of the 
system. To overcome this problem Berg and 
Neuhaus [2] have recently introduced multicanon- 
ical simulations which are based on reweighting 
ideas and can, in principle, be combined with any 
legitimate update algorithm. Using local update 
algorithms (Metropolis or heat-bath) it has been 
demonstrated in several applications [3] that the 
growth of autocorrelation times with system size 
is reduced to a power-law, r oc V a with a > 1. 
For the two-dimensional (2D) g-state Potts model 
values of a ss 1.3 have been reported for q = 7 [4] 
and q = 10 [5]. 

Since by construction the multicanonical en- 
ergy distribution is constant over the interesting 
energy range, a random walk argument would 
imply an exponent a = 1 for an optimally de- 
signed update algorithm. Here we present for 
Potts models a variant of the multicanonical ap- 
proach based on non-local cluster updates which 
turns out to be optimal in this sense [6]. Basi- 
cally the idea is to treat the cluster flips in the 
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first place and to reweight the bond degrees of 
freedom instead of the energy. 

2. ALGORITHM 

The basis of cluster algorithms [7] is the equi- 
valence Zpotts = Z^k = Zrc [8] of the spin repre- 
sentation of the Potts model (cr t =l, . . . , q) 



with the Fortuin-Kastelein (FK) 

Zfk = X] 5Z II [P<W^.i + S bij,o] 

{«i} {>>„} (ij) 

and random cluster (RC) representations 

z R c=Y,p B i Nci{b ' l}) -> B = zZ b ^ 

{bij} (ij) 
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where p = exp(/3) — 1. Here bi 3 =0 or 1 are bond 
variables and N c ({bij}) denotes the number of 
clusters. A cluster is a set of sites which are con- 
nected by active bonds 6^ = 1. By differentiating 
In Z with respect to (3 it is easy to see that the 
average of the energy E can be expressed in terms 
of the average of the number of active bonds B, 



d\nZ 
d(3 
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(B). 
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This suggests that the bond histogram Pb(B) 
should develop for (3 = (3 t ±6(3 & pronounced peak 
around B ^ = — ^jE 0: d, and in the vicinity of 



0.001 



0.000.") 




r,oo 



looo 



!.-)()() 



:,ooo 



.-,.-,oo 



oooo 



o.->oo 



\E\ 



Figure 1. Canonical energy and bond histograms 
for q = 7, L = 60, and (i = 1.292283. The bond 
histogram is plotted vs [(p + l)/p]B 



the transition point fit a double-peak structure 
similar to the energy histogram P e (E). In fact, 
as is illustrated in Fig. 1 for q = 7 and L = 60, a 
plot of Pb versus ^^-B is hardly distinguishable 
from P e versus \E\. For other values of q and L 
the comparison looks very similar. 

In terms of Pb the slowing down of canonical 
simulations is thus caused by the strongly sup- 
pressed configurations between the two peaks, 
analogous to the well-known argument for P e . 
To enhance these probabilities we therefore in- 
troduce in analogy to multicanonical simulations 
a "multibondic" partition function 

Z mu bo = ^2 X! LI b'WAij.i+'^j.oWC 8 )' ( 5 ) 

{"i}{bii}<ij) 

where Wb{B) = P b ~ 1 (B) between the two peaks 
and Wb(B) = 1 otherwise. Canonical expectation 
values can always be recovered exactly by apply- 
ing the inverse of the reweighting factor Wb(B). 
According to (5) we update the bi 3 as follows. If 
a t ^ dj then the bond bi 3 is never active and we 
always set b t j =0. If Oi = dj then we define 
B' = B — b^ and choose new values bi 3 = or 
1 with relative probability Wb(B') : ptVb(B'+ 1). 
Since for Wb(B) = 1 we just recover the original 
Swendsen-Wang algorithm we can now proceed in 
the usual way by identifying clusters of connected 
spins and choosing new random values 1 ... q in- 
dependently for each cluster. 



3. RESULTS 

To evaluate the performance of the multibondic 
(mubo) cluster algorithm we performed simula- 
tions of the 2D Potts model (1) with q = 7, 10 
and 20 on different lattice sizes for temperatures 
where the peaks of P e (E) have approximately the 
same height [6]. For comparison we ran stan- 
dard multicanonical (muca) simulations using the 
heat-bath update algorithm with the same pa- 
rameters. In each run we recorded 100 000 mea- 
surements of E and B in a time-series file. Be- 
tween the measurements we performed several 
lattice sweeps to ensure that the autocorrelation 
times in units of measurements and thus the ef- 
fective statistics of practically uncorrelated data 
was roughly the same in all simulations. 

For direct comparison with previous work we 
define flipping times 4t b ' p by counting the num- 
ber of update sweeps that are needed to travel 
from E < -E m in to E > E max and back. The cuts 
-Brain, max are chosen as the peak locations E ^(L) 
of P e (E) . In our simulations we have tested if E 
has passed the cuts after each sweep. Our results 
for t e ' p obtained in multicanonical and multi- 
bondic simulations for q = 7 and 10 are shown in 
the log-log plots of Fig. 2. Let us first concentrate 
on the results for q = 7 where we show for compar- 
ison also the data from previous multicanonical 
simulations [4] and from Rummukainen's hybrid- 
like two-step algorithm which combines micro- 
canonical cluster updates with a multicanonical 
demon refresh [9]. Both cluster update versions 
show qualitatively the same behavior and, for 
L > 20, perform much better than the standard 
multicanonical algorithm. From least-square fits 
to t e ' p = aV a we estimate am 1.3 for multicanon- 
ical heat-bath and a ss 1 for multibondic clus- 
ter simulations; see Table 1. Unfortunately, for 
q = 10 and 20 the situation is less favorable for 
the multibondic algorithm. While we still find an 



Table 1 

Results for the dynamical exponent a in multi- 
canonical and multibondic simulations. 



10 



20 



muca 
mubo 



1.27(2) 
0.92(2) 



1.32(2) 
1.05(1) 



1.26(1) 
1.09(1) 
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Figure 2. Log- log plots of autocorrelation times r E w of the energy vs lattice size for q = 7 and 10. 



exponent of a « 1, the prefactor a turns out to 
be so large that we can take advantage of this 
asymptotic improvement only for very large lat- 
tice sizes. As can be seen in Fig. 2, for g = 10 the 
cross-over happens around L = 50. Extrapolating 
to L = 100 we estimate that the multibondic al- 
gorithm would perform for this lattice size about 
1.5 times faster than the standard multicanoni- 
cal heat-bath. For q = 20 the same comparison 
clearly favors the standard algorithm for all rea- 
sonable lattice sizes - and we certainly cannot 
recommend the new algorithm for large q. 



4. CONCLUSIONS 

The multibondic cluster algorithm provides a 
combination of cluster update techniques with re- 
weighting in the random bond representation. It 
is technically not more involved than the multi- 
canonical approach and one lattice sweep takes 
about the same CPU time. This new algorithm 
is optimal in the sense that the exponent a is 
consistent with the optimal random walk value of 
a = 1 and it clearly outperforms the multicanoni- 
cal heat-bath algorithm. Compared to Rummu- 
kainen's intricate algorithm r B ip is smaller by a 
factor of ss 1.5. For larger values of q, however, 
the prefactor a turns out to be relatively large, 
rendering the new algorithm for reasonable lattice 
sizes more efficient than multicanonical simula- 
tions only for q<qo with go somewhat above 10. 
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